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Abstract 


This  report  contains  an  account  of  geostrophic  vort-^ics 
on  a  rotating  sphere.  A  vortex  of  this  type  3c  characterized 
by  a  singular  spherical  harmonic  of  degree  v  and  order  zero. 
This  function  has  a  role  analogous  to  that  of  the  Bessel 
function  which  characterizes  a  rectilinear  geostrophic  vortex 
in  a  rotating  plane. 


Preceding  page  blank 
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1.  Introduction 


A  thin  spherical  layer  of  an  incompressible,  inviscid  fluid 
which  is  held  on  the  surface  of  a  rotating  ball  by  gravitation 
can  be  taken  for  some  purposes  as  an  approximation  to  the  Earth's 
atmosphere.  The  analysis  of  the  two  dimensional  vortical  motion 
in  such  a  layer  should  be  useful  for  the  understanding  of  certain 
observed  meteorological  phenomena. 

In  a  previous  report  a  linearized  analysis  was  presented 
under  the  assumption  that  the  outer  surface  of  the  layer  is  a 
sphere  always  concentric  to  the  ball  representing  the  Earth. 

This  report  presents  a  linear  analysis  in  which  the  outer  surface 
of  the  fluid  layer  is  allowed  to  be  free  but  the  Coriolis  force 
is  assumed  to  be  independent  of  latitude.  This  assumption  de¬ 
fines  a  motion  which  can  be  regarded  as  a  first  approximation  to 
the  actual  motion  even  if  the  fluid  layer  is  not  confined  to  a 
narrow  band  bounded  by  two  circles  of  latitude. 

The  investigation  covers  the  case  in  which  tangential 
accelerations  are  neglected  in  comparison  with  the  Coriolis  force. 
In  other  words,  we  study  geostrophic  vortices  on  a  sphere  as 
contrasted  with  geostrophic  vortices  on  a  plane.  The  latter  have 
been  discussed  by  several  authors  in  connection  ^ith  a  tangential 
plane  approximation  to  the  motion  of  a  thin  layer  of  fluid 
covering  a  rotating  ball.  References  can  be  found  in  Morikawa, 
[1]. 
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2.  Equations  of  Motion 


Let  p  denote  the  distance  of  a  point  from  the  center  of  a 
ball  E  of  large  radius  a  which  rotates  with  constant  angular 
velocity  co  about  a  polar  axis.  Let  *  and  8  denote  respectively 
the  longitude  and  the  colatitude  of  a  point  on  the  rotating 
spherical  surface  S  of  E.  I*t  p  -  a  and  p  -  a  +  h(M,t)  repre¬ 
sent  two  surfaces  which  contain  an  incompressible,  Invlscid  fluid 
which  is  gravitationally  attracted  by  E.  Suppose  h  is  small 
compared  with  a. 

The  general  problem  is  to  find  the  motion  of  the  fluid 
after  its  constant  rotatory  motion  is  disturbed  by  the  sudden 

creation,  at  a  reference  time  t  =  0,  of  concentrated  vortices 
normal  to  S. 

The  velocity  of  a  fluid  particle  relative  to  S  is  defined 
by  the  components 

u  =  (p  sin  6)  ij|  =  tangential  component  toward  the  east; 
d© 

V  =  ~  P  cU  =  tangential  component  toward  the  north; 

w  =  =  radial  component. 

If  the  only  body  force  acting  Is  that  due  to  the  gravitational 

potential  0  of  E,  then  the  basic  hydrodynamical  equations  which 

define  the  motion  of  the  fluid  relative  to  S  are  the  continuity 
equation 

D  ^  +  pw  +  1  fdu  d(v  sin  6)7 

P  ^  *  ilnT  [5?  '  - L\  -  0 
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and  the  momentum  equations 


6Qp  sin  y  dij)  * 


In  these  equations  the  differential  operator  with  respect  to  the 
time  means 

d  _  d  .  u  d  _  v  £  ,  d 
cT£  "5T  p  sin  6  ^  p  "5^  "5p 

The  symbol  6Q  denotes  the  constant  density  of  the  fluid;  and  p 
stands  for  the  pressure  in 

9  p  p 

6  p  a)  sin  6 
Pl  *  p  +  60G - j - 

which  can  be  referred  to  as  the  modified  pressure. 

Since  h  is  supposed  to  be  small  compared  with  the  large 
radius  a  of  E;  and  since 

w((j>,e,a,t)  =  0 

let  us  neglect  the  radial  velocity  and  the  radial  variation  of  u 
and  v.  Let  us  assume  that  G  =  gp  where  g  is  a  constant;  and  that 
the  centrifugal  effects  manifested  by  the  partial  derivatives  of 

p  p  2 

6  p  o)  sin  6 

0  _  can  be  ignored.  Let  us  also  assume  that  the  motion 

is  such  that  the  nonlinear  terms  in  the  tangential  momentum 


~  -  u.v  -£££—£  +  —  +  2o)W  sin  6  -  2cdv  cos  6  =  - 
at  o  p 


dv  +  ™  +  u2  cot  e  +  2a)U  cos  e  _  1  |£i 

at  p  p  &0p 


12  ,  dp, 

«  -  - - - - 2a3U  sin  6  =  -  JL  -r-i  . 

at  P  P  60 
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equations  can  be  neglected;  and  that  the  radial  momentum  equa¬ 
tion  can  be  replaced  with  the  hydrostatic  law 

P(<M>p»t)  =  g6Q(h+a-p) 

which  satisfies  the  condition  that  the  pressure  is  zero  at  the 
free  surface  p  =  a+h. 

Under  the  above  assumptions  and  with  the  notation 
u(cj>,e,t)  =  u(  (|>,  £?,  a,  t)  , 
v(4>,  6,  t)  =  v(4>,  6,  a,  t )  ; 


an  approximation  to  the  motion  is  determined  by  the  equations 


(2.1) 


dw  _ 

"5p  a  s 


TnT  [’ 


du  d(v  sin  6) 


3- 


=  0  ; 


(2.2) 


-  2cdv  cos  6  =  - 


g  dh  . 
a  sin  0  ’ 


(2.3) 


If  +  cos  e  =  S  H  • 


An  integration  of  (2.1)  gives 


w(M,a+h,t)  -  a'sln  5 


d(  v  sin  6) 


W 


du 

"5^ 


The  kinematic  condition  at  the  free  surface  is 


w(  <|>,  6,  a+h,  t)  =  ^ 


and  hence 
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(2.4) 


dh  =  h 
cH;  a  sin  6 


5(v  sin  6) 


"50" 


.  & 
"5$ 


If  h(<t>6,0)  =  hQ  =  const.;  and  if  we  introduce 

h  -  h 

n(<M>  t)  =  — p —  > 


assuming  that  rj  is  small  almost  everywhere,  then  a  linearization 
of  (2.4)  yields 


(2.5) 


1  f5(v  sin  6)  5u 
^t  “  a  sin  0  50  c^i 


while  the  momentum  equations  become 

(2-6>  |u  -  2o£  COS  fl  -  -  ~  'g  ||  ; 

gh  n 

(2.7)  +  2cou  cos  6  =  — —  -g$  . 


3 .  Geostrophic  Vortices 

If  we  neglect  the  variation  of  the  Coriolis  force  with 
latitude  and  take 

03  cos  6  =  0)  cos  =  0)^ 


the  last  equations  of  Section  2  reduce  to 


(3.1) 


^t  a  sin  0 


5(v  sin  6)  5u 

50  ^ 
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(3.2) 

6u 

~SE  “ 

2o)jV 

(3.3) 

bv 

SE  + 

2a>^u 

gho  . 
a  sin  e  ^  ' 


gho  6 


a 


$ 


Hereafter  we  confine  ourselves  to  a  study  of  these  equations. 

As  will  be  explained  in  the  sequel,  they  lead  to  what  are  called 
geostrophic  vortices. 

The  elimination  of  q  from  (3.2)  and  (3.3)  leads  to 


(3.4) 

The  quantity 


d 

SE 


( 


1 

a  sin  0 


dn  sin  0  ,  6v 

— 55 —  +  5$ 


)  -  2m1ilt  . 


C  = 


a  sin  0 


d(u  sin  0)  ,  6v 

- 55 - + 


is  the  radial  component  of  vorticity;  and  by  integration  of  (3.4) 
we  have 


(3.5) 


1 

a  sin  0 


6(u  sin  0) 

- SB - 


20^1)  +C(^>,0,O)  . 


This  implies  that  we  can  use  (3.5),  (3.2)  and  (3.3)  as  a  basic 
set  of  equations  instead  of  (3.1);  (3.2)  and  (3.3). 

Equations  (3*2)  and  (3.3)  give 

ao,  |3  -  -j-a  ifa. 

1  60  sin  0  dtdcf) 

-  - 

and 


(3.6) 


2~ 
d  u 

it? 


,  2-  gh0 

4a)lu  =  — 
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(3.7) 


+  Jtafv  = 


'**1  a 

11  + 

a 

sin  0  15 

|  5tdS 

If  these  are  used  to  eliminate  u  and  v  from  (5«5  )>  we  find 


(3.8) 


sin 


(sln  6)  ^  + 


gH; 


^ +  U 


a  20)^  (<J),  6,  0) 
gfil 


For  a  vortex  of  strength  |i  concentrated  at  (<|>,,  ©1)  when 
t  =  0  we  take 


&22^(^e,0)  a22o)1  6(4)  -  4>1)6(©  -  6^ 

=  ^  a2sin  6. 


where  6  symbolizes  the  Dirac  delta  function. 

The  theory  of  Laplace  transforms  can  tu  used  to  show  that 
the  steady  state  solution  of  (3.8)  corresponding  to  a  con¬ 
centrated  vortex  is  such  that 


(3.9) 

where  f  must  satisfy 

(5*10)  sIFTT  W  (sin  e) 


2cd,^(<|>,  6 ) 
=  £ 


gn" 


1  S2*  ^a2oj2i/' 

TIFe  IF  "  “ 

M.6  ( 4>  -  <{>1)6(e  -  Si) 
sin  8. 


7 


The  function  if/  does  not  depend  on  t  and  the  associated  time 
independent  velocity  comonents  are 


(3.H) 

(3.12) 


U  “  Lt^oou^,  0,t)  “  f  ; 


V  -  Lfc^vt^O.t)  = 


1 

a  sin  S 


Note  that  (3.10)  is  a  consequence  of  the  equations  (3.11),  (3*12) 
and  the  vorticity  equation 

<3-13>  mriHTr  [|j  sln  ®)  +  U]  -  2“iLt-coT' 

m-6 (4>  -  4>1)5(e  -  ex) 

gho  a^sin  ©1 


In  other  words,  (3*10)  is  implied  by  (3.5),  (3.2)  and  (3*3)  when 
we  ignore  inertial  forces. 

The  motion  in  a  thin  planar  layer  of  fluid  tangential  to 
the  surface  of  the  Earth  is  often  used  as  an  approximation  to  the 
actual  motion  of  the  Earth's  atmosphere.  For  such  an  approxima¬ 
tion  the  analogues  of  (3.11);  (3*12)  and  (3.13)  are 


(3.14) 


(3.15) 

and 


(3.16) 


V 
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1  Sri'  ,  3V 

7  ST  + 


4m2x  u6($  -  iM6(r  -  r.) 

sV +  n 
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In  terms  of  the  polar  coordinates  (r,  <f>),  these  show  that  \  must 


satisfy 

(3.17) 


13  3*  ,  1  32*  n8(*-+1)8(r-r1) 

r-5?r3r  +  7^-  gV  =  - *1 


The  only  physically  admissible  solution  of  (3 .17)  is 


(3.18)  X  =  -  K0 


2  co,  /— p - p - 

- — -  -2rr1  cos  ( 4>  -  cf^ ) 


✓  Sh 


where  K  (  )  denotes  the  zeroth  order  modified  Bersel  function 

o 

of  the  second  kind.  This  defines  what  is  called  a  geostrophic 
vortex.  The  motion  of  various  configurations  of  such  vortices 
has  been  studied  by  various  authors,  notably  G.  K.  Morikawa  [1], 
[2]  whose  papers  contain  detailed  explanations  and  other 
references.  In  keeping  with  what  seems  to  be  accepted  termin¬ 
ology  we  can  say  that  (3.10)  defines  a  geostrophic  vortex  on  a 
sphere. 

Our  object  now  is  to  investigate  the  nature  of  the  solution 
of  (3.10).  This  equation  can  be  written  as 


<3,19)  si?TB  W  <sln  6) 


+  — V  — 4  +  v(v+u* 

sin  6  d<j> 


where 


(3.20) 


ii6(<j>  -if)1)6(e  -©]_) 
sin  0^ 


v (v+1) 


4co2a2 

^7  ’ 


The  solution  of  (3*19)  subject  to 


9 


a<e<£;  0  <  <t>  <  2tt  ; 

and  prescribed  boundary  conditions  on  the  circles  of  latitude 

0  =  a  i  9  =  p  ; 


namely  a  Green's  function  for  a  zone,  can  be  expressed  in  the 
form 


(3.21) 


[\i^(c°S  0)  +  Bm0j1(cos 


0 )  ]  cos  m(  (J>  -  <{)^ ) 


where  f^1  and  Q^1  are  solutions  of 


d 

cTz 


(1-  z2) 


v(v+l)  - 


m 


1  -z‘ 


f  =  0 


The  appropriate  series  (3.21)  can  be  summed  for  the  case  in  which 
the  zone  becomes  a  sphere.  However,  since  the  details  of  the 

analysis  based  on  the  use  of  (3.21)  are  lengthy,  let  us  proceed 
in  a  different  way. 

The  general  solution  of 


1  _  d 

sin  6  ‘5F 


(sin  0 ) 


v  (v+1  )i(/ 


0 


can  be  expressed  in  the  form 


V  =  c1pv(-cos  6)  +c2Pv(cos  9 ) 


when  v  is  not  an  integer  as  is  the  case  when  v  is  defined  by 


or 


v (v+1) 


4o)2a2 
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V 


where 


A  = 


In  the  neighborhood  of  the  north  pole  P^(cos  0)  is  continuous 
and  Pv(l)  =  1;  but 

P  (-cos  0) 
v v  ' 


behaves  like 

P  (-cos  0)  ou  2JlrP.  21..  i  n  0  . 

v  '  t r 


=  ?  ,  (-cos  0) 
-^  +  i  A 


00 


=  —  cosh  Air 
ir 


/ 

0 


cos  Axdx 


«/2(cosh  x  -  cos  0) 


This  suggests  that  if  a  vortex  is  concentrated  at  (fj^,  ©1);  if  d^ 
is  the  geodesic  distance  from  this  point  to  an  arbitrary  point 
(<J>,  0)  on  S;  and  if  is  the  angle  =  d-^/a,  then 


fr  -fe- -  P  (-cos  )  s  „ — £ -  P 

4  sin  vtt  v  1 '  4  sin  vtt  v 


M- 


-COS  0  COS  0- 


-sin  0  sin  0^  cos  (fjj-fj^) 


is  the  fundamental  solution  of  (j5.19)» 

It  can  now  be  verified  by  direct  substitution  that 


<3-22)  *  =  tts iH-nr  pv(-cos 

does  indeed  satisfy  (3*19)»  Furthermore,  computations  with  the 
velocity  components 
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“  4  a  sin  vtt  ^  V"C0S  7X)  i 


V  4a  sin  e  sin  vtt  Pv^"cos  7^)  > 

show  that  in  the  neighborhood  of  (<|).^,  t9^),  (3.22)  defines  a 
vortical  motion;  but  the  vortex  point  itself  remains  at  rest  —  it 
possesses  no  autonomous  motion.  In  fact,  it  can  be  shown  that 


Ja->GD  4" 


M- 

sin 


VTT 


Pv(-cos  7l) 


/ 


2  2 

r  +  r1  -  2r^ 


cos 


( 4>  -  4^) 


Hence  we  can  say  that  (3«22)  represents  a  geostrophic  vortex  on 
a  rotating  sphere  in  the  same  way  as  we  say  that 


X  “ 


K 


2  a). 


Ah; 


2rr1  cos  (<|>  - 


represents  a  geostrophic  vortex  on  a  rotating  plane. 

For  an  arbitrary  distribution  of  n  vortices  on  the  sphere 
the  function  if/,  which  can  be  regarded  as  a  stream  function,  is 


if  = 


M-i^y  (-cos 


where  7l  =  d^/a.  and  di  is  the  geodesic  distance  on  S  from  the 
point  of  concentration  of  the  i-th  vortex,  to  an 

arbitrary  point  ( <J>,  9).  That  is. 


cos  =  cos  9  cos  ei  +  sin  9  sin  6^  cos  (ij)  -  $  )  . 
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The  associated  velocity  field  is  given  by 


U  -  slf  -  4a  sin  vir  W  £=I  4iPv(-c°s  •y1)  ; 

V  “  a  sin  e  'jf  “  aa  sin  a  sin  vtt  dj,  4ipv  (-cos  yi )  . 

Since 

U  =  a  sin  6  ;  V  =  -a  || 


the  equations  of  motion  of  the  vortex  at  (<J>  ,Q  )  are 

K  K 


d(k 


"3T - 2 — ! - 

4a  sin  6  sin  vtt 


S  n 

'SU  11^  HiPv(-cos  7i) 

i?k 


and 


<M} 

6=0, 


d6k 


—  g 

4a  sin  0  sin  vtt 


n 


M'iPv(-C0S  V 

i^k 


H 

e=e. 


We  are  now  in  a  position  to  develop  a  general  theory  for 
geostrophic  vortices  on  a  sphere  analogous  to  the  basic  Helmholtz- 
Kirchoff  theory  for  rectilinear  vortices  normal  to  a  plane. 

Instead  of  doing  this  here  it  seems  more  useful  to  examine  a  few 
special  cases  which  may  be  approximately  applicable  to  certain 
physical  situations. 
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Geostrophlc  Vortices  in  the  Northern 


Hemisphere.  Either 


Velocity  Componant  Zero  at  the  Equator 

Let  us  turn  now  to  the  case  of  a  concentrated  geostrophlc 
vortex  in  the  northern  hemisphere  subject  to  the  condition  that 
the  normal  velocity  of  the  fluid  at  the  equator  is  zero.  For 
this  case  the  stream  function  ij/  must  satisfy 

sTFTT  W  (sin  +  v  (v+l )if/ 

sin  6  ^4> 

-  <t>1)6(e  -  e1) 

sin  y, 

for  1 

0  1  e,6l  <  ^  »  0  1  4*2.  ^  2tt  . 

The  boundary  condition  is 

(4,2)  V(*,§)  =  i  =  ° 

or 

2(1)^ 

glT  M'5)  =  ^  =  const. 

We  have  supposed  that  0)  =  0;  and  we  now  suppose  that 

nU^t)  =  0  for  all  t.  The  boundary  condition  (4.2)  then  be¬ 
comes 

(^•3)  ^(<j>,§)  =  0  . 

The  solution  of  (4.1)  which  satisfies  (4.3)  can  be  formed 
from  the  fundamental  solution 
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T-sln  vtt  pv(-cos  -Va) 


is  we  assume  that  the  idea  of  reflection  across  the  equator  will 
lead  toward  the  solution.  In  fact,  it  is  easy  to  verify  that 

f  pv  [-cos  0  cos  e±  -  sine  sin  Q±  cos  (<j)-<j)  )]': 

(4.4)  ^  =  ^ ; 

4  sin  vtt  )  > 

[-Pv  [cos  0  cos  ei  -  sine  sin  6±  cos  ($-<{>)] 


is  the  desired  solution.  With  respect  to  the  whole  sphere,  ip  of 
(4.4)  has  concentrated  vortices  at  (<^,0.,)  and  (4>1,tt-©  ) .  It 
should  be  noted  that  if  we  seek  a  stream  function  for  a  geo- 
strophic  vortex  on  a  rotating  plane  such  that  the  normal  velocity 
of  the  fluid  is  zero  along  a  circle  whose  center  is  at  the  center 
of  rotation,  then  a  result  analogous  to  (4.4  )  does  not  exist. 

The  velocity  components  of  the  motion  defined  by  (4.4)  are 


U. 

4a  sin  vtt 


f 

p' 

V 

-cos  0  COS  01 

• 

sin  0  cos  0-^ 

/ 

-sin  0  sin  0^  cos  ( 4> -  (jj^) 

-cos  0  sin  cos  (<J>  -t}^) 

+p’ 

V 

cos  0  cos  0^ 

• 

sin  0  cos  0^ 

V 

-sin  0  sin  0^  cos  ((j)-(j)^) 

+cos  0  sin  0^  cos  (ij>-  <j>^) 

and 
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-cos  0  cos  0 ^ 

0  sin  0 ^  cos  (<))  -  <J>^) 
cos  0  cos  0 ^ 

0  sin  0^  cos  ( <|>  -  (f>^) 


sin  0.^  sin  (<b  - 


sin  sin  ( 4>  -  ij^) 


J 


The  equations  for  the  motion  of  the  vortex  concentrated  at 
( (j^,  0^)  come  from  evaluating  at  ( e l  ^  velocity  field  due 
to  the  vortex  at  (<j>-pTr-0-^) .  They  are 


♦l- 


e,  = 


d^ 

cTT" 


d0x 

mr 


cos  6. 


2a  sin  vtt 


Py  (cos  2ex)  ; 


=  0 


From  these  we  see  that  the  vortex  remains  on  its  initial  circle 
of  latitude  which  it  traverses  with  constant  angular  velocity. 

The  function  ^  corresponding  to  a  concentrated  vortex  in 
the  northern  hemisphere  with  the  boundary  condition 


U(4>,§)  =  0 


can  also  be  found  by  appealing  to  a  reflection  process.  It  is 


*  ~  T'-sTE 


VTT 


P^t-cos  6  cos  ©1-  sin  0  sin  ^  cos  (<{>-  ) ] 

+Pv[cos  9  cos  9 ^  -  sin  0  sin  0^  cos  ( <|>  —  3 
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It  should  be  remarked,  however,  that  a  reflection  technique  for 
the  determination  of  solutions  of  (3.19)  for  a  general  zone  is 
unknown. 

If  the  northern  component  of  the  velocity  of  the  fluid  is 
zero  at  the  equator;  and  if  we  have  a  polar  vortex  of  strength  nQ 
with  three  others  of  equal  strength  p.  which  at  t  =  0  are  situated 
at 

(<J>p  ®i)  >  +  >  ( ^l  +  "5"*  ® l ) 

in  the  northern  hemisphere,  then  the  stream  function  is 


^  4  sin 


VTT 


+n 


O 


[ 

[ 

[ 


P  ( -cos  0  )  -  P  ( cos  © )] 

V  7  y  '  7  J 

-cos  0  cos  0.^  -  sin  0  sin  0 
cos  0  cos  0.^  -  sin  0  sin  0 


1 

1 


cos  (<j>-<|>^ 
cos  ( 4>— <f>1 


A  calculation  shows  that  the  polar  vortex  remains  stationary 
while  the  paths  of. the  other  three  are  determined  by 


(4.5) 


\  ‘  0 


f  ^0[PV(-C0S  )  +  Py  ( cos  ©! )] 


^k  “  -  / 

4a  sin  vrr  ' 


3  cos  0-^  (|  sin201  -  1) 
+(jL  -2  cos  ©.^(cos  20^) 

+  cos  0^(1  sin20]_) 


=  O  . 
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In  other  words,  the  position  of  the  vortices  subsequent  to  t  =  0 
are  given  by 

(<t>0>0)  ; 

(<|>1+ot,e1)  ;  (^1  +  ^L  +ot,e1)  ,*  (♦1+3t +ot,e1)  . 

Equation  (4.5)  shows  that  for  a  certain  strength  the  angular 
velocity  f)  becomes  zero  and  all  the  vortices  remain  stationary. 
The  stability  of  this  equilibrium  configuration  can  be  investiga¬ 
ted  by  using  the  methods  explained  in  the  paper  by  Morikawa  [1]; 
and  the  paper  by  Morikawa  and  Swenson  [2]. 
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